Topological vortex formation in a Bose-Einstein condensate 
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Vortices were imprinted in a Bose-Einstein condensate using topological phases. Sodium conden- 
sates held in a IofTe-Pritchard magnetic trap were transformed from a non-rotating state to one with 
quantized circulation by adiabatically inverting the magnetic bias field along the trap axis. Using 
surface wave spectroscopy, the axial angular momentum per particle of the vortex states was found 
to be consistent with 2% or 4/i, depending on the hyperfine state of the condensate. 
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As superfluids, Bose-Einstein condensates support ro- 
tational flow only through quantized vortices. The 
atomic velocity field is proportional to the gradient of 
the phase associated with the macroscopic wavefunc- 
tion. This phase winds through an integer multiple of 
2tt around a vortex line. Such a phase winding can be 
imprinted onto the condensate wavefunction either dy- 
namically or topologically. Dynamically, the phase of 
the condensate evolves according to the time integral of 
its energy, which can be tailored locally with a spatially 
varying external potential. Topologically, the phase of 
the condensate advances through adiabatic variations in 
the parameters of the Hamiltonian governing the sys- 
tem. This phase, which is solely a function of the path 
traversed by the system in the parameter space of the 
Hamiltonian, is known as Berry's phase 

In this Letter, we implement the proposal of Refs. 
||, [![ |[ and demonstrate the use of topological phases 
to imprint vortices in a gaseous Bose-Einstein conden- 
sate (BEC). Previously, vortices have been generated in 
two-component condensates using a dynamical phase- 
imprinting technique || and in single-component con- 
densates by rotating the cloud with an anisotropic po- 
tential |?], g, ^, [HJ, by slicing through the cloud with 
a perturbation above the critical velocity of the conden- 
sate [^, pj| , and through the decay of solitons [[l2| 
In this work, 23 Na condensates were prepared in either 
the lower, \F, ttif) = |1,— 1), or upper, |2,+2), hyper- 
fine state and confined in a IofTe-Pritchard magnetic trap. 
Vortices were created by adiabatically inverting the mag- 
netic bias field along the trap axis and could be removed 
by returning the bias field to its original direction. Us- 
ing surface wave spectroscopy ]l4|, 15, [l6|, we measured 
the axial angular momentum per particle of the |1, — 1) 
and |2, +2) vortex states to be consistent with —2mph 
as predicted [||||. 

A IofTe-Pritchard magnetic trap consists of an ax- 
ial bias field (with curvature) and a two-dimensional 
quadrupole field in the orthogonal plane [[l7[ |l8[ |l9| : 

B(x, y, z) = B z z + B\xx - yy), (1) 
where B' is the radial magnetic field gradient and 



quadratic terms have been neglected. For a condensate 
of radial extent R, inverting B z from B z ^> B'R > 
to B z <C —B'R < rotates the atomic angular momen- 
tum, F, by 7r radians. While all atomic angular momenta 
rotate through the same angle, a relative phase is es- 
tablished across the condensate because the angular mo- 
menta rotate about a unit vector h(<p) = sin <\> x + cos <p y 
that depends on the azimuthal angle, </>, describing the 
atomic position (Fig. 0(a)). 

As B z is inverted, F adiabatically follows B(x, y, z) 
and the condensate always remains in the state \F, mp) 
with respect to the local magnetic field. However, in a 
basis fixed in the lab frame, the condensate makes the 
transition \F,m z = -fra^) — > \F,m z = — ttif), where 
ttif and m z are the projection of F along the local mag- 
netic field direction and z-axis respectively. Applying the 
quantum mechanical rotation operator gives the conden- 
sate wavefunction in the lab frame after inverting B z as 
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where T is the angular momentum operator such that 
F — (T) and p(f) is the number density of condensed 
atoms. The topological phase factor e - l2rn F4> describes a 
vortex of winding number 2|mi?| with the sense of rota- 
tion dependent on the sign of mp. 

This result can be interpreted in terms of Berry's 
phase Fig. |l|(b) shows the orientation of F in the lab 
frame for mp > 0. Atoms located at position k = 1,2 
in Fig. [|(a) have azimuthal angle and angular mo- 
mentum Fk. As B z is inverted, Fk traces path k from 
top to bottom on the sphere in Fig. 0(b). The topolog- 
ical phase acquired by an atom in this process is solely 
a function of the path traversed by its angular momen- 
tum vector. Since this path depends on the azimuthal 
angle describing the atomic position, a relative phase is 
established between spatially separated atoms. The con- 
densate wavefunction after inverting B z is given by 
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FIG. 1: Geometry of the rotating magnetic field for the 
topological vortex formation. (a) The unit vectors 6(0) 
(solid arrows) point in the direction of the two-dimensional 
quadrupole field providing the radial confinement of a Ioffe- 
Prit chard magnetic trap. The atomic angular momenta rotate 
about the unit vectors n(0) (open arrows) as the axial bias 
field, B z , is ramped from positive to negative values, (b) For 
an atom in state \F, ttif), its atomic angular momentum, F, 
traverses a path on a sphere of radius as it adiabati- 

cally follows its local magnetic field. The primed coordinate 
system is centered on the atomic position and has axes par- 
allel to those of the unprimed coordinate system in (a). For 
an atomic position described by the azimuthal angle 0, F 
rotates in a half-plane defined by 0' — —<j> for tuf > and 
0' = — + 7T for vrtF < as B z is inverted. After inverting B z , 
the relative topological phase acquired between atoms located 
at positions 1 and 2 in (a) is proportional to the solid angle 
subtended by the shaded surface, bounded by the contour 
marked with arrowheads (see text). 



where 7(0) is the topological phase acquired by atoms 
with azimuthal angle 0. 

For an atom in state |F, m^), Berry's phase, 7(C), ac- 
quired as its angular momentum vector traverses a closed 
contour, C, on the surface of the sphere in Fig. |l|(b) is 
given by 

7 (C) = -m F Sl(C), (5) 

where Q(C) is the solid angle subtended by a sur- 
face bounded by the contour C. The relative phase, 
7(^1 ) — 7(^2), is unaltered by closing the contours traced 
by each Fk along an arbitrary but identical path. For 
clarity, we choose to close each contour along path 2 it- 
self and hence 7(0i) — 7(^2) = 7(C), where C is the 
contour formed by path 1 traversed from top to bottom 
and path 2 traversed from bottom to top, as indicated 
with arrowheads in Fig. |l|(b). 

A surface bounded by this contour subtends a solid 
angle Q(C) = 2(0i — 02), yielding a relative phase 7(0i) — 
7(^2) = — 2771^(01 — 02). Thus, we make the assignment 

7 (0) = -2m F 0, (6) 

up to an additive term independent of position. This 
is the same phase as in Eq. [| reinterpreted in terms of 
Berry's phase. 



In this work, Bose-Einstein condensates containing 
over 10 T 23 Na atoms were created in the |1,— 1) state 
in a magnetic trap, captured in the focus of an optical 
tweezers laser beam, and transferred into an auxiliary 
"science" chamber as described in Ref. ]20| . While opti- 
cally confined by the tweezers, condensates were prepared 
in the |2,+2) state by sweeping through the |1,— 1) <-> 
|1,0) <-> |1,+1) radio-frequency transition with 100% ef- 
ficiency, then sweeping through the |2,+2) 
microwave transition with 80% efficiency In the 

science chamber, the condensate was loaded into a mi- 
crofabricated IofTe-Pritchard magnetic trap formed by a 
Z-shaped wire carrying current I and an external mag- 
netic bias field, B± : as detailed in Ref. [p2| . Condensates 
were detected via axial absorption imaging whereby reso- 
nant laser light propagating along the z-axis illuminated 
the atoms and was imaged onto a CCD camera [[l9| . 

Typical wiretrap parameters were / = 1200 mA, 
B± = 5.4 G, and B z w 1 G, resulting in a radial mag- 
netic field gradient of B' = 120 G/cm. For atoms in 
the |1,— 1) state, the axial and radial trap frequencies 
were u z = 2ir x 6 Hz and co± = 2ir x 210 Hz respec- 
tively. For atoms in the |2,+2) state, the wiretrap fre- 
quencies (for identical magnetic field parameters) were 
larger by a factor of \[2 due to the larger magnetic mo- 
ment. After transfer into the wiretrap, condensates in the 
|1, -1)(|2, +2)) state contained over 2 x 10 6 atoms(l x 10 6 
atoms) and had a lifetime in excess of 10 s(3 s) with an 
applied radio- frequency shield Jl9| ]. This represents the 
first magnetic trapping of 23 Na condensates in the upper 
hyperfine level, with previous work done exclusively in 
optical dipole traps [ pl[ . 

Along the wiretrap axis, the magnetic field is 

B(x = 0,y = 0,z) = (B z + \b"z 2 )z, (7) 

where quadratic terms neglected in Eq. p] have been in- 
cluded. The axial magnetic field curvature, B" which 
arises from the geometry of the Z-wire, was held con- 
stant throughout the experiment. By reversing the ex- 
ternal axial magnetic field, we inverted B z . Changing the 
sign of B Z1 but not B" resulted in a magnetic field saddle 
point at the center of the cloud and axial anti-trapping 
of weak- field seeking atoms. This limited the condensate 
lifetime after inverting B z to < 50 ms. 

Vortices created by inverting B z were identified by 
characteristic centrifugal density depletions observed af- 
ter ballistic expansion (Fig. ^|). These vortices could be 
removed by returning B z to its original direction. 

For condensates in the |1,— 1) state, the best results 
were achieved by inverting the axial bias field linearly 
from B z = 860 mG to —630 mG in 11 ms. For atoms in 
the 1 2, +2) state, the optimum ramp time over the same 
range was 4 ms. The field inversion process caused an 
atom loss of « 50% due to non-adiabatic spin-flips as B z 
passed through zero, in reasonable agreement with cal- 
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FIG. 2: Observation of vortices formed by imprinting topo- 
logical phases. Axial absorption images of condensates in the 
|1, —1) state after 18 ms of ballistic expansion (a) prior to 
inverting B z , after inverting B z and holding the trapped con- 
densate for (b) 5 ms and (c) 20 ms, and (d) after inverting 
B z and then returning it to its original direction. Axial ab- 
sorption images of condensates in the |2, +2) state after 7 ms 
of ballistic expansion (e) prior to inverting B z , after invert- 
ing B z and holding the trapped condensate for (f) ms and 
(g) 5 ms, and (h) after inverting B z and then returning it to 
its original direction. The field-of-view is (a)-(d) 570 /xm x 
570 fim and (e)-(h) 285 /im x 285 /im. 



dilations in Refs. ||, |||. The density depletions shown 
in Fig. ||(b,c,g) were observed after inverting the ax- 
ial bias field and holding the trapped condensate for 
longer than a radial trap period. Thus, the atom loss 
from the center of the cloud during the field inversion 
process could not be responsible for the observed den- 
sity depletions. Typical |1, — 1)(|2, +2)) condensates af- 
ter inverting the axial bias field contained up to 1 x 10 6 
atoms(0.5 x 10 6 atoms) with a Thomas-Fermi radius 
Rtf — 5.4 ± 0.2 iivh^Rtf = 4.5 ± 0.2 /im) in a trap with 
radial frequency co_\_ = 2tt x 250 Hz(cu± — 2ir x 350 Hz). 

The axial angular momentum per particle of the 
vortex states was measured using surface wave spec- 
troscopy |l4|, [l5[ [l6| . A superposition of counter-rotating 
(rri£ = ±2) quadrupolar (£ = 2) surface waves was ex- 
cited in the condensate by radially displacing the mag- 
netic trap center for 200 /is. Here, t and rri£ charac- 
terize the angular momentum and its projection along 
the z-axis of the quadrupole modes respectively. This 
created an elliptical condensate cross-section with time- 
dependent eccentricity. In the absence of a vortex, the 
vfi£ — ±2 quadrupole modes are degenerate and the axes 
of the elliptical condensate cross-section remain fixed in 
time. This degeneracy is lifted by the presence of a vor- 
tex, causing the axes to precess in the direction of the 
fluid flow. The precession rate, 6, is given by |l4|, |l5[ [l(J 



6 = 



(L z ) 
2M(r 2 _) 



(8) 



where (L z ) is the axial angular momentum per particle 
characterizing the vortex state, M is the atomic mass, 



By measuring the precession rate of the quadrupole 
axes and the mean-squared radius of the condensate in 
the trap, the axial angular momentum per particle was 
determined. After exciting the quadrupolar modes, the 
condensate evolved in the trap for variable times in the 
range 0.2 — 7.4 ms. The condensate was then released 
from the trap and imaged with resonant light after bal- 
listic expansion as shown in Fig. ||(a)-(l). The result- 
ing images were fit to an elliptical Thomas-Fermi pro- 
file to determine the orientation of the quadrupole axes. 
The orientation angle is plotted as a function of time in 
Fig. ||(m). To determine the mean-squared trapped con- 
densate radius, vortices were imprinted in the conden- 
sate but quadrupolar modes were not excited. Images 
of ballistically expanded condensates similar to those in 
Fig. |||(b,c,f,g) were fit to a Thomas- Fermi profile with cir- 
cular cross-section. The fitting routine ignored the cen- 
tral region of the cloud where the density was depleted 
due to the vortex core. The mean-squared trapped con- 
densate radius was derived through the relation 
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and (rj_) 



2 \ 



y 2 ) is the mean-squared trapped con- 



densate radius with vortices present. 



where R± is the Thomas- Fermi radius of the condensate 
after ballistically expanding for a time r from a trap with 
radial frequency uo±. The factor 1 + uo^r 2 accounts for 
the change in Thomas-Fermi radius during the expansion 
process [^3| and the factor 2/7 results from averaging 
over the inhomogeneous condensate density distribution 
assuming no vortices are present. For low angular mo- 
mentum vortex states, the density depletion at the vortex 
core does not significantly modify the density distribution 
of the condensate and we expect the 2/7 factor to still 
be accurate Q. 

For condensates in the |1,— 1) state, the quadrupole 
oscillation was excited after a delay of 0, 5, and 20 ms 
from the completion of the inversion of the axial bias 
field. The measured axial angular momenta per parti- 
cle were +1.9(2)(2)ft, +2.1(2)(2)ft, and +1.9(l)(2)ft re- 
spectively, where the first uncertainty is associated with 
the linear fit to the precession angle and the second un- 
certainty is associated with the determination of (r*j_). 
For condensates in the |2,+2) state, the quadrupole os- 
cillation was excited immediately upon the completion 
of the inversion of the axial bias field. The measured 
axial angular momentum per particle was — 4.4(1) (4)7i. 
For both internal states, the measurements are consistent 
with the predicted axial angular momentum per particle 
of-2m F ft|,|g|. 

Multiply charged vortices are unstable against decay 
into singly charged vortices (25). From our experiments, 
we cannot determine if the condensate contained one 
multiply charged vortex or multiple, singly charged vor- 
tices. If multiple vortices were present, they must be 
closely spaced since they were not resolved after ballistic 
expansion. Furthermore, if the singly charged vortices 
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FIG. 3: Surface wave spectroscopy. Axial absorption images 
after 18 ms of ballistic expansion of 1 1 , — 1) condensates under- 
going a quadrupole oscillation (a)-(d) in the presence of a vor- 
tex and (e)-(h) in the absence of a vortex. Successive images 
were taken during successive half periods of the quadrupole 
oscillation such that the short and long axes of the elliptical 
cross-section were exchanged. Images (a)-(d) show counter- 
clockwise (positive) precession of the quadrupole axes, while 
images (e)-(h) show no precession, (i)-(l) Axial absorption 
images after 7 ms of ballistic expansion of |2, +2) condensates 
undergoing a quadrupole oscillation in the presence of a vor- 
tex. The images were taken during a single half period of the 
quadrupole oscillation. Images (i)-(l) show clockwise (neg- 
ative) precession of the quadrupole axes. The field-of-view 
is (a)-(h) 570 /xm x 570 /xm and (i)-(l) 285 /xm x 285 /xm. 
(m) Precession angle vs. time in the presence of a vortex for 
|1, —1) condensates measured after a delay of ms (open cir- 
cles), 5 ms (open squares), and 20 ms (open triangles) from 
the completion of the inversion of the axial bias field, in the 
absence of a vortex for |1, —1) (open diamonds) and |2, +2) 
(filled diamonds) condensates, and in the presence of a vortex 
for |2, +2) condensates measured immediately upon the com- 
pletion of the inversion of the axial bias field (filled circles). 



had moved apart considerably, it would have lowered the 
extracted value of (L z ) [|l5|, which was not observed even 
with delayed probing. 

In conclusion, we have used topological phases to im- 
print vortices in a Bose-Einstein condensate. The salient 
feature of this phase imprinting technique is the passage 



of the zero of a magnetic quadrupole field through the 
condensate. Confining the condensate optically would 
allow for a series of such passages and would lead to 
states with very high angular momentum. This tech- 
nique opens the potential for studying the stability of 
multiply charged vortices and the dynamics of vortex- 
vortex interactions at short separations. 
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